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Abstract 

In a number of classical statistical-physical models, there exists a characteristic dimensionality 
called the upper critical dimension above which one observes the mean-field critical behavior. In- 
stead of constructing high-dimensional lattices, however, one can also consider infinite-dimensional 
structures, and the question is whether this mean-field character extends to quantum-mechanical 
cases as well. We therefore investigate the transverse-field quantum Ising model on the globally 
coupled network and on the Watts-Strogatz small-world network by means of quantum Monte 
Carlo simulations and the finite-size scaling analysis. We confirm that both of the structures ex- 
hibit critical behavior consistent with the mean-field description. In particular, we show that the 
existing cumulant method has difficulty in estimating the correct dynamic critical exponent and 
suggest that an order parameter based on the quantum-mechanical expectation value can be a 
practically useful numerical observable to determine critical behavior when there is no well-defined 
dimensionality. 

PACS numbers: 05.30.Rt,75.10.Jni,89.75.Hc 



1 



The effects of geometry on physical properties have drawn attention from various areas 
of statistical physics. First of all, we note the extensive studies on complex networks (see 
Ref. |l for a review). Aside from their ubiquity and practical importance, a remarkable point 
is that many of the complex networks tend to enhance correlations dramatically when a 
statistical-physical system is put on top of them. Let us look at the Watts-Strogatz (WS) 
network,^ for example: it starts from a one-dimensional (ID) ring of size where each point 
is linked to its 2k nearest neighbors. Among the total Nk bonds, we pick up PNk bonds with 
< P < 1, and then for each of them, we change one of its endpoints to a randomly chosen 
site. The WS network is a classical model of a small-world network characterized by two 
structural features: the length scale L is only logarithmic with respect to the system size A^ 
and the clustering coefficient is nevertheless relatively high.- Due to the small-world property, 
one observes the mean-field critical behavior in many cases, such as percolation, Ising, and 
XY spin models.-"- Since that is expected only for very-high-dimensional structures, many 
complex networks, including the WS network, are often called infinite dimensional. One may 
notice that this terminology also makes sense in that L ~ N^/^ in dimensional lattices 
while the WS network has only L ~ logA^. The theory of finite-size scaling has been well 
developed for such complex networks in combination with extensive numerical calculations.- 

Compared to the recent progress in the study of classical systems in these structures, 
however, there is relatively little known about quantum cases. Even the transverse-field 
quantum Ising system on the Bethe lattice has been examined very recently, and it is still 
to be clarified whether the phase transition is really of the mean-field type,-"^ and the same 
model on the WS network has been checked only in terms of the thermal phase transitions.— 
In this work, we thus investigate how the quantum Ising spin system in the transverse field 
behaves in the infinite-dimensional globally coupled network and the WS network. We do 
not consider the scale-free network in this work since the critical behavior can be possibly 
far richer than the conventional mean-field description.— The transverse-field Ising model 
is defined by the quantum-mechanical Hamiltonian 

where ai = (af , af , erf) is a spin-^ operator at each given site i and the summation runs over 
all the nearest-neighbor pairs. This Hamiltonian also contains a constant J > to represent 
the ferromagnetic-interaction strength and another constant F to represent the transverse- 



field strength. If located in d dimensions, this system exhibits a phase transition at a certain 
critical-field strength, F = Fc, and its criticality is known to be equivalent to that of the 
classical {d + l)-dimensional Ising model where F is replaced by temperature T.— But it 
is not entirely obvious what is going to happen when the dimensionality lacks its precise 
meaning in the infinite-dimensional structures, such as the WS network. Before studying 
the WS network, however, it would be good to check another infinite-dimensional structure, 
i.e., the globally connected network, since the transverse-field Ising system on this structure 
has been already known to exhibit the zero-temperature quantum phase transition of the 
mean-field universality by the Hamiltonian-matrix diagonalization.— The globally coupled 
network is simply obtained by connecting every site to all of the others, but the Hamiltonian 
needs a little modification to Eq. ([T]) in order to make it extensive as follows: 

H = -N-'jY^ata^^+TY^at. (2) 

i^j i 

The main tool in analyzing Eqs. ([1]) and ([2]) is the world-line quantum Monte Carlo method,— 
where the system is endowed with an additional temporal dimension in the imaginary-time 
direction. The length in this direction is the same as the inverse temperature h = 1/T 
and divided into Nt slices. The Suzuki- Trotter decomposition says that we get correct 
results of the original transverse-field Ising model at temperature T a.s Nt — ?■ oo. Our 
quantum Monte Carlo algorithm is based on the single-cluster flip algorithnt^ and employs 
the idea of continuous imaginary time, as described in Ref. |l9|. For convenience, we actually 
work with integer variables in programming (instead of floating variables in indexing) the 
imaginary-time axis, and then make the range of the indices, A^^, up to O(IO^). This virtually 
eliminates the error involved in the width of an imaginary-time slice in the Suzuki- Trotter 
decomposition. We have checked that the implemented algorithm convincingly reproduces 
the expected results for d = 1 and 2 where the dynamic critical exponent is already known 
to be z = 1.— The dynamic critical exponent z tells us how to scale the temporal size b with 
respect to L to observe correct critical behavior, so it means that 6 ~ = in ordinary 
d-dimensional lattices. Note that this quantity is required since our quantum Monte Carlo 
algorithm cannot directly access T = 0, but needs a suitable extrapolation of T in relation 
to the system size, unless one simulates extremely low T compared to all of the finite sizes 
under consideration. When z is known a priori, everything is just straightforward. Applying 
the algorithm to the infinite-dimensional structures described above, however, we usually 



do not know it at the starting point and need to estimate z' to scale T ~ 'N~^' in a suitable 
way. We can formulate the problem as follows: let us consider the magnetic order parameter 
\m\ taken from the whole [d + l)-dimensional space-time region, that is 



{|m|> = iVj-'/V-' 



E 



i,t 



(3) 



where each time slice is indexed by t and the bracket (■ ■ •) means the thermal average, as 
well as the disorder average in the case of the WS network. Note that a-^^ in Eq. ([3]) is not 
an operator, but a spin value projected onto the z direction at a certain space-time. One 
can argue a finite-size scaling ansatz in the following two-parameter form: 



m\) = N-^/^'f \{T - T,)N^/''' ,TN'' 



(4) 



It is usually complicated to directly deal with this two-parameter scaling form. But when 
T is suitably scaled with a dynamic critical exponent z' so that r = TN^' = const., Eq. (jl]) 
reduces to a usual scaling form with a single parameter, from which it is possible to determine 
the critical exponents. We therefore follow the procedure proposed in Refs 
estimate z'. The idea is that Binder's cumulant for a finite-sized system. 
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and 
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to 



U = 1 



(5) 



3 (m^)^ ' 

vanishes not only in the limit of T — )■ oo, but also in the limit of T — )■ where the system 
virtually becomes a classical ID Ising chain elongated in the imaginary-time direction. One 
may thus expect that a maximum will appear by f/ = f/max at a certain characteristic 
temperature Tmax for each given system size A^. If the field strength is above the critical 
threshold by F > Fc, then the peak value f/max will eventually vanish as becomes larger, 
but will approach a nonzero value if F < Fc instead. This speculation provides a way to 
determine F,, where f/max remains constant regardless of A^. We illustrate how it looks for 
the ID case at Vi^^^ = 1 in Fig. [I](a). Note the nice symmetry as the size and inverse 
temperature b vary at the same time. It implies that one can also assume that U obeys a 
finite-size scaling form as 



U = g 



(F-FjA^i/^',TA^^' 



(6) 



Therefore, f/ is a function of r = TN^' at F = F^, and a correct z' will yield a scaling 
collapse of U [Fig. [Dl^b)]. It is also notable that f/max in Fig. [1] is fully consistent with 
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FIG. 1: (Color online) (a) Binder's cumulant for the ID transverse-field Ising model at F = T^^^ = 
1 as a function of the inverse temperature b = 1/T. The transverse- field strength T is in units of J 
and the temperature T is in units of J/ks where ks is the Boltzmann constant, (b) Scaling collapse 
with the dynamic critical exponent z = 1. The dotted lines indicate its universal amplitude ratio 
for the two-dimensional Ising model [U* = 0.610692(2)]. 

the universal amplitude ratio U* = 0.610692(2) for the two-dimensional Ising model. By 
checking the size dependence of 6max = 1/^max where U reaches the maximum, one will 
be able to obtain the dynamic critical exponent. Then, we calculate the magnetic order 
parameter with changing both N and T, and finally get the critical exponents in the zero- 
temperature thermodynamic limit by the standard finite-size scaling, ~ N~^/'^' and 
dU/dV ~ N^l"' at F = F^. 

In addition, it is interesting to note that in Fig. [U^b) the curve exhibits refiection symme- 
try with the peak position as the axis of symmetry. The reason for this can be argued in the 
following way: for given L, let us begin with observing the system at 6 = fomax- Then one 
can say that the correlation length in the spatial direction amounts to L, while that in the 
temporal direction amounts to h. Let us increase the temperature so that h becomes 6max/2. 
When a system has different lengths in different axes, the correlation length will be bounded 
by the shortest one, so the correlation length will become about one half compared to the 
previous case of 6 = femax- It means that as we traverse the system in the spatial direction, 
we will find no correlation at a distance around -L/2, and we may roughly regard the system 
as separated into two uncorrelated parts in the spatial direction. If we instead lower the 
temperature to h = 26max, the correlation length still remains bounded by L, and we will 
find no correlation as we traverse the system by 6max in the temporal direction. In this 
sense, we may regard the system as separated into two uncorrelated parts in the temporal 
direction this time. It becomes precisely symmetric to the case of 6 = &max/2 by exchanging 
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the axes and adjusting the overall scale. In other words, there is symmetry between the 
cases of logfemax — ^ log&max ± log 2, which explains the reflection symmetry in Fig. [IJb). If 
we extend this argument to a higher dimension, say d > 1, the system will be broken into 
2°^ pieces when we take b — )■ &max/2. Since the other part of the argument about b — ?■ 26rnax 
remains the same, it predicts that U will decay faster on the low-6 side than on the high-6 
side, and that this asymmetry will be more pronounced as d gets higher. 

Let us check how the method using Eq. ([6]) works for the globally coupled network 
with Eq. ([2]). As mentioned above, it has been already established that the transverse- 
field quantum Ising model on it exhibits the mean-field critical behavior, and the critical 
threshold is at = 1 in units of J.— Since the transverse-field Ising model has the upper 
critical dimension (i„ = 3, it also leads to z' = z/d^ = 1/3.— We thus calculate Binder's 
cumulant at F = 1 and depict the results in Fig. [2]^a). We immediately notice two features 
from the figure: first, the maximum height is close to the universal amplitude ratio for the 
mean-field model, U* ~ 0.270520.— Second, the curve shapes become more asymmetric as 

gets larger. This observation is qualitatively explained by the argument in the paragraph 
above, and furthermore implies that it is not possible to make a scaling collapse of this plot 
just by tuning T in accordance with as postulated in Eq. ([6]), at least up to the sizes that 
we have in this work [Fig. [2](b)]. It hinders one from using the extrapolation to A^ — > oo 
and T — 7- from the finite-sized data at finite temperatures by Eq. ([6]). Even in the case 
that one works at a fixed temperature, the finite-size effects may well enter and dominate 
the situation: when the temperature is relatively high, say T = 0.5, we have little difficulty 
finding the mean-field transition at Fc ~ 0.958(5) with the c/asszca/ upper critical dimension 
^(classic) _ value of U is close to the universal amplitude ratio [Fig. [2](c) and its 

scaling collapse (not shown here)]. This observation can be explained as follows: rewriting 
Eq. as U = g{TN' ) at F = F^, one may expect g to converge to a constant when 
the system gets large enough, i.e., A^ ^ T^^/^'. It in turn means that the scaling form in 
Eq. becomes insensitive to the second argument at a sufficiently high T. This is why we 
find such finite-size scaling behavior with varying F at T = 0.5 in Fig. [2](c). The classical 
mean-field behavior at this T clearly confirms a consequence of crossover phenomena from 
the quantum to the classical critical behavior at different temperatures.-^ However, if we 
move to a low but still finite temperature, i.e., T = 0.1 in this example, the finite-size effect 
becomes so substantial that it gets hard to characterize the transition [Fig. [2](d)]. Note that 
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the crossing point is still moving to the low-F side, rising up very slowly to U*. If we increase 
by more than an order of magnitude, Binder's cumulant begins to move in the correct 
direction [Fig. [2](e)]. The critical-field strength at T = 0.1 can be estimated to differ from 
r = 1 by less than 5 x 10~^, and the cumulant will behave as t/ = f/*(l + cN~^ + ■ ■ ■) 
with a coefficient c and a correction exponent A > at this critical point. Since Eq. 
leads to U = g{bN^^') at F = Fc, one may guess that A = z' = 1/3 by expanding the 
scaling function g into a Taylor series. This simple guess is not very far from our numerical 
observation shown in Fig. [2](f). To sum up, our results suggest that Eq. ([5]) might not be a 
very efficient numerical observable to reproduce the correct results for infinite-dimensional 
structures. 

As discussed above, although intuitively appealing, the concept of the characteristic tem- 
perature Tmax where U shows a maximum is not adequate for high- dimensional structures. 
The difficulty can be overcome by using an alternative quantity to U. Let us take the 
magnitude of spontaneous magnetization at an arbitrary fixed point in the imaginary-time 
axis. 



{\s\) = N-' 



(7) 



Due to the time translational symmetry, defined in Eq. ([7]) is a time-independent value, 
and, moreover, it corresponds to the Boltzmann statistics with eigenf unctions diagonalizing 
the Hamiltonian Eq. ([T]) or Eq. ([2]), i.e., 

I '12 

where Z is the partition function, and \ sj\ = \aj\ \{si \ \si) \ and Ej are the quantum- 

mechanical expectation value of the operator ^ • erf and the energy eigenvalue for the jth 
eigenstate = J^i^i respectively. Consequently, as T — )■ 0, goes to a nonzero 
expectation value of the ground state for finite N, which is different from the behavior of 
(|m|) in Eq. ([3]). The scaling form is again considered as 

= N~^/''h \{T - T,)N^/''', TN''] , (8) 

so one finds scaling collapse {\s\) N^^'^' = h{TN^') at F = F^. Furthermore, as mentioned 
above. Fig. [3](a) shows that this quantity actually converges to a constant at 6 — i- oo for each 
A^. Since is insensitive to T as T — )■ 0, we can drop out the second argument of the 
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FIG. 2: (Color online) Observations for the globally coupled network, (a) Binder's cumulant at 



r = r 



(G) 



1. (b) Scaling of the horizontal axis with the known dynamic critical exponent z' = 1/3 



[compare with Fig. mb)]. (c) Binder's cumulant at T = 0.5 shows a crossing point, already close 
to the universal value, while (d) the crossing is still moving due to finite-size effects at T = 0.1. (e) 
Binder's cumulant as a function of at three different F's, with temperature T = 0.1 as in (d). 
The dotted lines represent the universal amplitude ratio for the mean-field case [U* ~ 0.270520]. 
(f) The difference of U from U* at T = 0.9999, in comparison with N~^^^. 

scaling function h in Eq. ([8]) at sufficiently low T and we are back to the single-parameter 
scaling form, 

{\s\) ^ N-^^"' h \{T -T,)N^/'''' 

around Tc- Fitting the values of a.t T ^ 1.4 x 10~^ as a function of N, for example, 
gives us an estimate of fS/u' = 0.339(2), which is fairly close to 1/3. Then we can estimate 
z' ^ 0.3 by making the data for different A^'s collapse into a single curve at large b [see, e.g., 
Fig.[3](a)]. One might say that the mean-field behavior could have been found by using 
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FIG. 3: (Color online) Observations for the globally coupled network, (a) Magnetization at a 
fixed imaginary time [Eq. ([7])] yields scaling collapse with z' = 1/3. (b) Scaling of the magnetic 
order parameter [Eq. ([3])] with z' = 1/3 at T = 1. By scaling T with this z' in such a way that 
r = TN^' = 2, (c) Binder's cumulant and (d) the magnetic order parameter show the correct 
mean- field quantum phase transition at Ec =1 with j3/i>' = 1/3 and 1/u' = 2/3. The dotted line 
represents the universal amplitude ratio for the mean-field case \U* « 0.270520]. 

instead [see Fig. [3](b)]. But our point is that is easier to handle in the sense that it 
converges a well-defined value in the Iarge-6 limit, which allows us to get around the problem 
of working with two scaling parameters at the same time. Beginning with scaling (|s|), one 
finds numerical observables such as U and {\Tn\) consistent with the mean-field behavior 
[Figs. |3](c) and [31(d)] . Therefore, we suggest that Eq. ^ can be a practically useful guidance 
in measuring the dynamic critical exponent. The advantage of using f/, on the other hand, 
comes from its universal values containing a good deal of information on universality classes, 
and it indeed works well in estimating Fc. Figures |3]^a) and [3]^b) furthermore suggest that 
a relevant scaling region should be found at 6 ^ 1, rather than around the peak position. 
Once again, we note that U vanishes at this large 6, while does not vanish for finite A^, 
so that the latter quantity can be more easily handled at the scaling region. 

Let us now proceed to the WS network. We choose = 3 and P = 0.05 and take dis- 
order averages over independent realizations. The number of realizations varies but usually 
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FIG. 4: (Color online) Observations for the WS network, (a)-(c) Binder's cumulant as a function 
of b at three values of T. In each panel, the curves represent = 200, 400, 800, 1600, and 3200 
from left to right, respectively, (d) From Eq. ([6]), we measure the dynamic critical exponent 
z' = 0.197(8), whereby this scaling collapse is obtained. The dotted lines represent the universal 
amplitude ratio for the mean-field case [U* ~ 0.270520]. 



amounts to O(IO^). We begin with the existing method using U, assuming that it will work 
here as in ID. Then we find characteristic temperatures for various sizes of the WS network 
by measuring U, from which we estimate ri^^"* = 4.72(1). [Figs. Hl^a) toHl^c)]. Noting that 
the small sizes tend to overestimate U, we determine F^ by observing the largest sizes. It 
is interesting that we do not see such large asymmetry as in the globally coupled network 
[compare Fig. |2]^a) with Fig. Hl^d)], which might be due to the fact that the length scale of 
the WS network is a logarithmic function of N. Another interesting indication is that the 
height of ?7max at F = 4.72 is close to the mean- field universal amplitude ratio [Figs. Hl^b) 
andlH^d)]. The dynamic critical exponent at F = 4.72 is estimated as z' = 0.197(8) from 
the dependence of T^ax on [Fig. Hl^d)]. It is apparently very close to 1/5, and leads to 
f3/u' = 0.296(8) and = 0.598(16), or, in other words, (3 = 0.49(3) and z/' = 1.67(4). 
Although they are close to the mean-field values 1/2 and 3/2, respectively, we have no 
explanation for z' = 1/5. 

Since we already know that Eq. ([6]) does not always give correct z', we may instead try 
the alternative method relying upon Eq. (|8]). By fitting the values of (|s|) at T ^ 1.9 x 10^^ 
with different A^'s, we find P/u' = 0.33(1), which contains 1/3 within the error bar. The 
dynamic critical exponent is then estimated as z' ^ 0.3 by making the data collapse [see, 
e.g.. Fig. El^a)]. This method definitely prefers z' = 1/3 to 1/5, which implies that the use 
of Eq. ([6]) might yield a spurious estimate of z' as in the globally coupled network. As 
seen in Fig. |5]^a), our observations of fS/u' and z' are fully consistent since both of them are 
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FIG. 5: (Color online) Observations for the WS network, (a) Scaling collapse of (|,s|) at F = 
-p(WS) _ ^ ^2 according to Eq. 1^, with p/v' = 1/3 and z' = 1/3. (b) Scaling collapse of the 
magnetic order parameter with z' = 1/3. Choosing this z' , we find (c) the crossing of U and (d) 
scaling collapse of at F = 4.72 with the mean-field values, i.e., /3/z^' = 1/3 and l/i^' = 2/3. 

The dotted line indicates the universal amplitude ratio for the mean-field case [U* « 0.270520]. 

explained by = 3. Note also the nice agreement in Fig. Mh)- Accepting this z' = 1/3, 
the crossing point of U is also consistent with the universal amplitude ratio of the mean- 
field type, and the best fits with and dU/dV give us /3 = 0.50(3) and v' = 1.38(4) 
[Fig. EI^c)]. The value of u' is reasonably close to 3/2 and our numerical results do not rule 
out this value [Fig. EJ^d)]. Overall, we may consider our Monte Carlo results for the WS 
network as supporting the mean-field behavior, i.e., /3 = 1/2, u' = 3/2, and U* ~ 0.270520. 
It is worth noting that Eq. ([H]) again outperforms Eq. in the finite-size scaling analysis. 

In summary, we have studied the transverse-field Ising model on infinite-dimensional 
structures by quantum Monte Carlo simulations. Since the world-line Monte Carlo algorithm 
cannot directly access T = 0, we need the dynamic critical exponent z' for each of the 
structures. For the WS network, we have estimated z' = 0.3(1) and fS/u' = 0.33(1) from 
or P = 0.50(3) and z/' = 1.38(4) from and dU/dV under the assumption that 

z' = 1/?). These strongly support the mean- field behavior with the upper critical dimension 
du = 3, i.e., (3 = 1/2 and v' = 3/2. Drawing this conclusion, we have also found that the 
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existing finite-size scaling method using U can give incorrect values of z' in the infinite- 
dimensional structures. For this reason, we have suggested an alternative observable 
based on the quantum-mechanical expectation value, which is found to converge more quickly 
to the correct scaling behavior than U does. 

Acknowledgments 

We thank Hangmo Yi for discussions and comments. This work was supported by a Na- 
tional Research Foundation of Korea (NRF) grant funded by the Korea government (MEST) 
(Grant No. 2011-0015731). This research was conducted using the resources of the High Per- 
formance Computing Center North (HPC2N). We thank R. Oriis and J. Vidal for informing 
us of relevant works on the Lipkin-Meshkov-Glick model.— 



* Corresponding author. E-mail: |beomjun@skku.edu| 

^ S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Rev. Mod. Phys. 80, 1275 (2008). 

2 D. J. Watts and S. H. Strogatz, Nature (London) 393, 440 (1998). 

3 C. Moore and M. E. J. Newman, Phys. Rev. E 62, 7059 (2000). 

^ B. J. Kim, H. Hong, P. Holme, G. S. Jeon, P. Minnhagen, and M. Y. Choi, Phys. Rev. E 64, 
056135 (2001). 

5 H. Hong, B. J. Kim, and M. Y. Choi, Phys. Rev. E 66, 018101 (2002). 

^ K. Medvedyeva, P. Holme, P. Minnhagen, and B. J. Kim, Phys. Rev. E 67, 036118 (2003). 

H. Hong, M. Ha, and H. Park, Phys. Rev. Lett. 98, 258701 (2007). 

^ D. Nagaj, E. Farhi, J. Goldstone, P. Shor, and I. Sylvester, Phys. Rev. B 77, 214431 (2008). 

9 F. Krzakala, A. Rosso, G. Semerjian, and F. Zamponi, Phys. Rev. B 78, 134428 (2008). 

° H. Yi and M.-S. Choi, Phys. Rev. E 67, 056125 (2003). 

1 H. Yi, Eur. Phys. J. B 61, 89 (2008). 

2 H. Yi, Phys. Rev. E 81, 012103 (2010). 

3 R. J. Elhott, P. Pfeuty, and C. Wood, Phys. Rev. Lett. 25, 443 (1970). 
^ P. Pfeuty, J. Phys. C 9, 3993 (1976). 

^ J. Um, S.-I. Lee, and B. J. Kim, J. Korean Phys. Soc. 50, 285 (2007). 



12 



R. Botet and R. JuUien, Phys. Rev. B 28, 3955 (1983). 

N. Kawashima and K. Harada, J. Phys. Soc. Jpn 73, 1379 (2004). 

U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

H. Ricger and N. Kawashima, Eur. Phys. J. B 9, 233 (1999). 

S. Sachdev, Quantum phase transitions (Cambridge University Press, Cambridge, 1999). 

M. Quo, R. N. Bhatt, and D. A. Huse, Phys. Rev. Lett. 72, 4137 (1994). 

H. Rieger and A. P. Young, Phys. Rev. Lett. 72, 4141 (1994). 

J. Salas and A. D. Sokal, J. Stat. Phys. 98, 551 (2000). 

E. Luijten and H. W. J. Blote, Int. J. Mod. Phys. C 6, 359 (1995). 

W. A. C. Erkelens, L. P. Regnault, J. Rossat-Mignod, J. E. Moore, R. A. Butera, and L. J. 
de Jongh, EPL 1, 37 (1986). 

S. Dusuel and J. Vidal, Phys. Rev. Lett. 93, 237204 (2004); P. Ribeiro, J. Vidal, and R. Mosseri, 
ibid. 99, 050402 (2007); R. Oriis, S. Dusuel, and J. Vidal, ibid. 101, 025701 (2008). 



13 



